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Abstract. A new scenario is presented for the produc- 
tion of 23 Na in asymptotic giant branch (AGB) stars. The 
scenario takes advantage of the periodic third dredge-up 
episodes characterizing those stars, which mix primary 12 C 
from their intershell layers to their surface. Two successive 
intcrpulse/pulse/dredge-up sequences are then required to 
produce 23 Na. During the first sequence carbon and oxy- 
gen are converted into 14 N by the hydrogen burning shell, 
and subsequently transformed into 22 Ne by the helium 
burning shell. During the second sequence, 22 Ne is con- 
verted into 23 Na by the hydrogen burning shell, which 
is brought to the surface by the subsequent dredge-up 
episode. The 23 Na produced by this scenario is thus pri- 
mary. 

The efficiency of this scenario is analyzed through stan- 
dard evolutionary AGB model predictions combined with 
synthetic calculations for the surface chemical evolution. 
It is shown that primary 23 Na can efficiently be produced 
as soon as the surface C+N+O abundance enhancement 
reaches a certain level depending on the stellar metallicity. 
The required surface C+N+O abundance enhancement 
amounts to ~0.4 dex in solar metallicity stars, and to 
~0.8 dex at a metallicity five times less than solar. 

An analytical study of Na production further reveals 
that the surface 23 Na abundance asymptotically evolves 
to a 'line of primary sodium enrichment' (LOPSE) in the 
[C+N+O] -[ 23 Na] diagram. That LOPSE represents the 
23 Na abundance evolution predicted in zero metallicity 
AGB stars experiencing third dredge-up episodes. An an- 
alytical relation for the surface 23 Na abundance evolution 
as a function of the surface C+N+O abundance is pro- 
vided. 

The predicted surface 23 Na enhancements can exceed 
0.5 dex depending on the level of surface 12 C enrich- 
ment, and increases with decreasing stellar metallicity. 
The quantitative prediction of 23 Na surface abundances, 
however, is presently subject to a high level of uncertainty, 
partly due to the still poor quantitative prediction of the 
structural evolution of AGB stars (dredge-up episodes in 
particular) , and partly due to the uncertainties still affect- 
ing some nuclear reaction rates (such as 23 Na destruction 
by proton capture). 



The case of massive AGB stars in which hot bottom 
burning occurs is also discussed. The production of sec- 
ondary sodium in those stars is a natural consequence of 
22 Ne burning in their envelope, if the temperature at the 
base of the envelope exceeds 70 million K. It requires, how- 
ever, many interpulscs to be significant. The production 
of primary sodium from the dredge-up of primary 22 Ne 
and its subsequent burning in the envelope, on the other 
hand, is estimated not to be very efficient, expect maybe 
in low-metallicity stars. 

An eventual detection of high Na overabundances in 
carbon stars or related objects would support the sce- 
nario of primary sodium production in AGB stars. Such 
an observational evidence may have been found in at least 
one post- AGB star. Further observations of those objects 
are called for. Observations of 23 Na in planetary nebulae 
are also encouraged. Finally, the production of primary 
23 Na by AGB stars, if confirmed observationally, may have 
played a non-negligible role in the chemical evolution of 
our Galaxy. 

Key words: stars: AGB and post- AGB - stars: carbon - 
stars: evolution - stars: interiors - stars: abundances 



1. Introduction 

Sodium overabundances have been observed in A-F su- 
pergiants, with 0.0 <J [Na/H] <J 0.5 for the F supergiants 
([X] = logX star - logX ) and 0.7 £ [Na/H] £ 0.8 for the 
A-type supergiants (Takeda & Takada-Hidai 1994 and ref- 
erences therein). The presence of enhanced sodium abun- 
dances at the surface of those stars is explained by the op- 
eration of the Ne-Na mode of H-burning during the main 
sequence phase, which leads to an increase of 23 Na in the 
deep stellar layers by a factor up to ten. The penetra- 
tion in those layers of the convective envelope when the 
star becomes a red giant, process called 'first dredge-up', 
then brings the synthesized sodium to the surface. First 
dredge-up predictions indeed confirm the production of 
23 Na to the level observed in A-F supergiants (El Eid & 
Champagne 1995 ). 
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High 23 Na overabundances have also been measured 
at the surface of red giants in globular clusters as early as 



in the late 70's (see Mowlavi 1998 for references). Again, 
the Ne-Na chain of H-burning and the operation of first 
dredge-up are made responsible for those surface over- 
abundances. The observed overabundances, however, are 
higher than those predicted by standard/] model predic- 
tions, but can be explained if extra-mixing is assumed to 
operate below the convective envelope, probably due to ro- 
tationally induced meridional circulation. In both globular 
clusters' red giants and near solar metallicity supergiants, 
the production of 23 Na results from the transformation of 
the initial 22 Ne. The sodium produced in those stars is 
thus secondary^. 

In this paper, a new scenario is proposed to produce 
primary sodium in asymptotic giant branch (AGB) stars. 
The scenario, explained in more details in Sect. 0, takes 
advantage of the 12 C and 22 Ne production in the helium 
burning shell (HcBS) of those stars, and of the periodic 
mixing of those nuclei to the stellar surface. Part of the 
22 Ne may then subsequently transform into 23 Na in the 
hydrogen burning shell (HBS) of those stars. Such a sce- 
nario is expected to operate in carbon stars. 

Several standard model stars are computed from their 
pre-main sequence up to the AGB phase. The stellar 
masses and metallicities covered in those calculations 
range from M = 1.5 to 6M Q and from Z = 0.004 to 
0.02, respectively. The properties of those standard mod- 
els are presented in Sect. g. In particular, the surface 23 Na 
abundance evolution prior to the AGB phase, and the 
abundances of 12 C, 22 Ne and 23 Na predicted in the inter- 
shell layers during the AGB phase are discussed. Primary 
sodium production is then studied in detail in Sect. ||| us- 
ing synthetic AGB calculations. The case of the 6 M star, 
which undergoes hot bottom burning (HBB), is considered 
separately in Sect. ||. Finally, the observational implica- 
tions/confirmations of the scenario are briefly discussed 
in Sect. ||. Conclusions are drawn in Sect. 0. 

Two Appendixes are further presented at the end of 
the paper. Appendix A describes briefly the CNO and Ne- 
Na modes of hydrogen burning, and presents the yields of 
those elements at hydrogen exhaustion as a function of 
the H-burning temperature. In particular, the impact of 
the new NACRE reaction rates on the 23 Na production 
is discussed. Appendix B, on the other hand, presents an 
analytical study of 23 Na production in AGB stars. 



1 Standard models refer, in this paper, to models using the 
Schwarzschild criterion to delimit convective zones without ap- 
plying any extra-mixing procedure such as overshooting or dif- 
fusive mixing induced by rotation 

2 An element is said to be secondary if its synthesis requires 
the presence of some elements heavier than helium in the ini- 
tial stellar composition. In contrast, an element is said to be 
primary if it can be synthesized in a star of the first generation 
with zero metallicity (the metallicity being the mass fraction 
of all elements heavier than helium) . 
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Fig. 1. Structural evolution of the intershell regions of the 
2.5 M model star with Z = 0.008 during the 15th inter- 
pulse and 16th pulse. The origin of the abscissa is arbi- 
trarily set to the time of onset of pulse 16. Filled regions 
correspond to convective zones. The dotted lines identify 
the location of maximum energy production in the H- 
burning (top) and He-burning (bottom) layers. The differ- 
ent phases of an instability cycle and zones A, B and D in 
the intershell layers described in Sect. @ of the main text 
are also indicated. The simulation of a third dredge-up, 
as predicted by models using the prescription of Mowlavi 
( 1993 ) , is shown in the afterpulse phase (region hatched at 
45 degrees in the right subpanel). Note that the time ab- 
scissa has different scales in the two subpanels displaying 
intcrpulse 15 and pulse 16. 



2. AGB stars and 23 Na production 

2.1. AGB stars 

An AGB star is characterized by an electron degenerate 
C-0 core of 0.5 to 1.2 M , by a thin HcBS of about 
1O~ 2 M capped by a thin HBS of a few 10~ 4 M Q , and 
by a deep convective envelope extending from above the 
HBS up to the surface. The HeBS is known from numer- 
ical simulations to be thermally unstable and to liberate, 
periodically and on a short time-scale (of several tens of 
years), 10 2 to 10 6 times the energy provided by the H- 
burning shell. These outbursts lead to the development of 
a pulse convective tongue in the HeBS. Figure GL for ex- 
ample, displays the structural evolution of the intershell 
layers of a 2.5 M , Z = 0.018 model star between its 15th 
and 16th pulse. 

A thermal instability cycle is divided, for our purposes, 
into three phases (illustrated in Fig. m: the interpulse, the 
pulse and the afterpulse phases. The main characteristics 
of each of these phases in relation with sodium production 
are described as follow (a more detailed description can 
be found in Mowlavi & Meynet, in preparation, hereafter 
called MM99): 



Nami Mowlavi: Sodium production in asymptotic giant branch stars 



Interpulse phase: During the interpulse phase, hydrogen 
burns in a thin layer located at mass fractions increasing 
with time, while the HcBS is almost extinct. The ashes left 
over by that HBS accumulate in the underlying He-rich 
zone (called zone B, see Fig. pi). Among them, 14 N is the 
main product emerging from the CNO mode of hydrogen 
burning, while 23 Na is produced in the NeNa chain (see 
Appendix A). 



Pulse phase: As the mass of the He-rich zone increases, a 
thermal instability is triggered in the HeBS, and a pulse 
convective tongue develops. The pulse starts at the lo- 
cation of maximum energy production in the HeBS, and 
extends outward up to close the HBS. As it grows, it en- 
gulfs the ashes of H-burning available in zone B. 14 N, in 
particular, is mixed down to the high temperatures char- 
acterizing the pulse, and burns through cv-capture into 
18 0. This nuclide, in turn, may burn and synthesize 22 Ne 
through 18 (a , 7) 22 Ne . Those products are left over by 
the pulse in zone D, together with primary 12 C synthesized 
by the 3-cv reaction. 



Afterpulse phase: During the afterpulsc phase, the struc- 
tural evolution of the intershell layers is dominated by 
thermal relaxations. The temperature in those intershell 
layers drop, and the convective envelope deepens into the 
H-depleted and, eventually into the C-rich layers. The 
ashes of H-burning which escape the ingestion by the 
pulse convective tongue (zone A in Fig. gj), and part of 
the material of zone D processed by the pulse (containing 
in particular 12 C, 22 Ne and 23 Na) are mixed to the surface. 
This process of mixing C-rich material into the envelope 
is called 'third dredge- up' (3DUP). 



2.2. Sodium production 

Sodium is produced in the HBS through 22 Ne (p , 7) 23 Na . 
This reaction essentially transforms into 23 Na all the 22 Ne 
supplied by the envelope to the HBS. For solar initial 
composition (displayed in Table ||), it would lead to an 
intershell 23 Na mass fraction of 1.61 x 10~ 4 . At tempera- 
tures exceeding 35 x 10 6 K, the slight burning of the more 
abundant 20 Ne nuclei contributes to an extra-production 
of 23 Na (maximum 60% at 55-60 x 10 6 K, see Fig. || of 
Appendix A). Above 60 x 10 6 K, however, 23 Na may be 
destroyed by p-capture through 23 Na (p , a) 20 Ne and/or 
23 Na (p , 7) 24 Mg . A more detailed description of the NeNa 
chain is provided in Appendix A. The actual predictions 
from standard AGB model calculations are discussed in 



Table 1. Solar abundances, in mass fraction, of the most 
abundant isotopes of the C, N and O elements, and of 
the stable nuclei involved in the Ne-Na chain (the solar 
metallicity is taken equal to Z = 0.018). 



Sect. 3.2 



nuclei 


solar mass fraction 


12 C 


2.878 x 10 -3 


14 N 


1.049 x 1(T 3 


16 


9.103 x 10~ 3 


20 Ne 


1.536 x 10" 3 


21 Ne 


3.925 x 10~ 6 


22 Ne 


1.236 x 10 -4 


23 Na 


3.166 x 10~ 5 



It must be noted that the impact of that intershell sec- 
ondary 23 Na on the surface sodium abundance is negligible, 
though. This is due to the fact that the mass of intershell 
material dredged-up to the surface is small (< 0.01 M©, 



see Fig. [l]) compared to the mass contained in the enve- 
lope U1M ). 

If we consider, however, the 12 C and 22 Ne dredged- 
up from the HeBS into the envelope, and the resulting 
nucleosynthesis during the interpulse and pulse phases, 
then a much larger amount of 23 Na may be synthesized. 
The chain of events is summarized as follow: 

1) 3DUP mixes into the envelope some primary 12 C pro- 
duced in the intershell layers by the 3-ev reaction; 

2) during the following interpulse, the 12 C (and 16 0) pro- 
vided by the envelope to the HBS is burned into 14 N 
through the CNO cycle, which accumulates in zone B of 
the intershell layers; 

3) the 14 N left over in zone B is engulfed in the next 
pulse, and contributes to the nucleosynthesis occurring 
in that pulse. In particular, 22 Ne is synthesized through 

14 N (a , 7) 18 F (/3+ ) 18 (a , 7) 22 Ne ; 

4) part of that 22 Ne emerging from the pulse is dredged-up 
in the envelope during the afterpulse phase; 

5) the 22 Ne provided by the envelope to the HBS during 
the next interpulse is transformed into 23 Na by the NeNa 
chain. The fraction of that 23 Na which survives destruc- 
tion by p-capture accumulates in zone B of the intershell 
layers; 

6) the 23 Na left over in zone B is engulfed in the next pulse, 
but is unaffected by He-burning nucleosynthesis. It thus 
emerges from the pulse unburned but diluted over zone D, 
and is partly mixed to the surface by the following 3DUP. 

In summary, the sodium produced in this chain of 
events essentially results from the conversion of the enve- 
lope's CNO into 22 Ne during an interpulse/pulse/3DUP 
chain of events, followed by the conversion of 22 Ne into 
23 Na during a second interpulse/pulse/3DUP chain of 
events. Two dilution processes are thus necessary. Assum- 
ing a solar repartition of the elements heavier than helium, 
it can be shown (Sect. IO) that the initial CNO abun- 
dances are insufficient to produce 23 Na in excess of its 
initial abundance. The dredge-up of primary 12 C from the 
HcBS is thus essential to increase the envelope's C+N+O 
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Table 2. 23 Na surface mass fraction predicted after first 
dredge-up {X{) and at the end of the E-AGB phase (A 2 ). 
The initial mass fraction Xi n u is also given for compar- 
ison. The mass fractions are given in units of 10 -5 . The 
increase of the surface 23 Na at the onset of the first pulse 
relative to its initial abundance is given in the fifth (in 
number ratio) and sixth (in dex) columns. 



model star 


-A-init 


Xi 


x 2 


X2 / Xinit 


M/M Q ,Z 


10 -5 


10~ 5 


10~ 5 


ratio 


dex 


6.0, 0.020 


3.52 


5.38 


6.32 


1.80 


0.25 


4.0, 0.018 


3.17 


4.82 


4.88 


1.54 


0.19 


3.0, 0.020 


3.52 


5.20 


5.20 


1.48 


0.17 


2.5, 0.018 


3.17 


4.54 


4.54 


1.43 


0.16 


1.5, 0.018 


3.17 


3.37 


3.37 


1.06 


0.03 


2.5, 0.008 


1.41 


2.17 


2.17 


1.54 


0.19 


2.5, 0.004 


0.70 


1.15 


1.16 


1.66 


0.22 



abundance, and to lead to significant 23 Na production. 
The sodium produced in this way is thus primary. 

3. Standard models 

3.1. The models 

Seven model stars are evolved from the pre-main sequence 
phase up to the 15-30th pulse in the AGB phase. The 
initial masses M and metallicities Z are (M j M Q , Z) = 
(1.5, 0.018), (2.5, 0.018), (3, 0.020), (4, 0.018), (6, 0.020), 
(2.5, 0.008) and (2.5, 0.004). 

The stellar evolution code originates from Mowlavi 
( |1995| , see also MM99 for more details about the mod- 
els). It was designed to follow as accurately as possible the 
structural and chemical evolution of AGB stars. The abun- 
dance profiles of 47 nuclei from H to S, linked by a nuclear 
reaction network of 160 reactions, are followed consistently 
with the stellar structure all along the stars evolution. The 
initial mass fractions of hydrogen and helium are taken to 
be X = 1 - Y - Z and Y = 0.23 + 2.5 x Z, respectively 
The abundances of the other elements are taken from An- 
ders & Grevesse ( |1989| ) scaled to the required metallicity. 
In particular, the solar abundances, taken at Z = 0.018, 
of those nuclei involved in the synthesis of 23 Na are listed 
in Table [j] The rates of the nuclear reactions involved in 
the Ne-Na chain are given in Appendix A. Recently, a new 
compilation of reaction rates has been made available by 
the NACRE project. Their impact on the production of 
23 Na in AGB stars is also discussed in that Appendix. 

The convective zones are delimited by the standard 
Schwarzschild criterion, with no extra-mixing beyond the 
convective borders so defined. Such standard models have 
been shown to fail to consistently reproduce the 3DUP 
phenomenon (see Mowlavi 1999] and references therein). 
The use of some sort of extra-mixing is required to model 
the 3DUP episodes. Such procedures, however, are very 
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Fig. 2. a Temperature and b abundance profiles of the 
stable Ne-Na nuclei, as labeled on the curves, in the inter- 
shell layers at the end of the 6th interpulse of the 2.5 M Q , 
Z = 0.018 standard model star. The extensions of zones A 
and B described in Sect. and of pulses before and after 
the considered interpulse are also shown in the lower part 
of panel b. 



computer time consuming. For that reason, and in order to 
analyze the sodium production in a variety of AGB models 
of different masses and metallicities, the models presented 
in this section are performed without extra-mixing. The 
surface 23 Na abundances are then computed a posteriori 
using synthetic AGB calculations. This is done in Sect, m 

3.2. Surface Na abundances prior to the thermally 
pulsing AGB phase 

In standard models not experiencing 3DUP episodes dur- 
ing the AGB phase, the surface 23 Na abundance of red gi- 
ants can only be affected by the first and, possibly, second^ 
dredge-ups. The surface sodium abundances predicted af- 
ter those dredge-up episodes are shown in Table for all 
our stars. The 23 Na abundance enhancements reach 40 to 
80% (0.15-0.25 dex) relative to its initial abundance in all 
the intermediate-mass stars considered in this paper, and 
less than 10 % in the 1.5 M Q star. Those predictions agree 
with the ones available in the literature. 



3.3. Intershell 23 Na abundances during the AGB phase 

As recalled in the preceding section, no surface 23 Na abun- 
dance alteration is predicted by our standard TP-AGB 



3 The second dredge-up occurs at the end of the early AGB 
phase (i.e. before the thermally pulsing AGB phase) and only 
for stars more massive than 
Z = 0.004). 
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Fig. 3. Same as Fig. g, but at the end of the 22nd inter- 
pulse of the 4M Q , Z — 0.018 standard model star. 



modelsr]. It is, however, instructive to analyze the abun- 
dance profile of the secondary 23 Na in the intershell layers 
of the standard models. The question of 23 Na destruction 
by p-capture, in particular, is of relevance for our study. 
The intershell temperature and Ne-Na abundance pro- 
files at the end of the 6th interpulse of the 2.5 M©, 
Z = 0.018 model star are shown in Fig. @. Secondary 
23 Na is produced at the level expected from the operation 



of the NeNa chain (Sect. 2.2). The temperatures charac- 
terizing the HBS are lower than 50 x 10 6 K, and no 23 Na 
destruction by p-capture is obtained, as expected. The in- 
tershell 23 Na abundance is 5 times higher than its surface 
abundance. 

The same profiles are presented in Fig. ||, but for the 
22nd interpulse of the 4M , Z = 0.018 model star. The 
HBS of those models is characterized by temperatures 
reaching 75 x 10 6 K, leading to the activation of 23 Na 
burning by p-capture. Indeed, a slight destruction of 23 Na 
is observed in the ashes left over by the HBS (Fig. ||). How- 
ever, the intershell 23 Na abundance is found to still be 2.3 
times higher than its surface abundance. Compared to the 
intershell abundance in the 6th interpulse of the 2.5 M©, 
Z = 0.018 star, the 23 Na destruction amounts to a fac- 
tor of about 2. This Na destruction factor agrees with the 
yields of the parametric H-burning calculations presented 
in Appendix A. 

Recently, a new compilation of nuclear reaction rates 
has been provided by the NACRE project (Arnould et 
al. 1999Q , which predicts much higher 23 Na p-capture rates 
than estimated before. In order to evaluate the impact of 
those new reaction rates on Na production in AGB stars, 




Fig. 4. 12 C mass fraction emerging from the pulses of our 
model stars as a function of the H-depleted core mass. 
Each mark locates a pulse. 



the 22nd interpulse of the 4M©, Z = 0.018 is recomputed 
with those new rates. The resulting intershell 23 Na abun- 
dance amounts now to only 70% of its surface abundance. 
The destruction factor is thus 7, which is much higher than 
with the old reaction rates. It must be noted, however, 
that the high uncertainties still affecting the p-capture 
rates do not, at present time, enable to lead to a definite 



conclusion. Arnould et al. (1999) show, for example, that 



4 The eventual dredge-up of the He-rich (but not C-rich) 
material from zone A (Fig. FEl) leads to a negligible alteration 
of the surface Na abundance. 



23 Na can still be produced to an abundance 3 times higher 
than its initial abundance, even at 80 x 10 6 K, within the 
uncertainties of the new nuclear reaction rates. 

3.4- Intershell 12 C and 22 Ne abundances 

The mass fraction of 12 C emerging from the pulses is 
shown in Fig. |]for all our model stars as a function of the 
H-depleted core mass. It reaches up to 30% in the 2.5 M©, 
Z = 0.018 star, and evolves after a sufficient number of 
pulses to an asymptotic value of ~ 20% independent of 
the stellar mass and metallicity. The abundance of 22 Ne 
emerging from the pulses, on the other hand, is shown in 
Fig. |k as a function of the temperature of the HBS at the 
onset of each pulse. 

Let us consider the 22 Ne production in the intershell 
layers. If we assume that all the available 12 C and 16 in 
the HBS are transformed into 14 N, then 1.24 x 10~ 2 mass 
fractions of 14 N are expected in the intershell layers for so- 
lar initial chemical composition (cf. Table [j]). Converted 
into 22 Ne, this would lead to a 22 Ne abundance emerging 
from each pulse, Xo( 22 Ne), 157 times higher than its ini- 
tial abundance X, n , t ( 22 Ne). For solar initial composition, 
this corresponds to Xo( 22 Ne) = 0.019 mass fractions. 

In reality, 12 C and 16 are present at the level of few 
percents in the ashes of H-burning, reducing accordingly 
the amount of 14 N available for the synthesis of 22 Ne. 
Though the reduction is very small (few percents), it is in- 
structive to pursue further the discussion. Parametric nu- 
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Fig. 5. a Abundance of 22 Ne, relative to its initial stellar 
abundance, emerging from the fifth pulse on in our model 
stars as a function of the temperature of the H-burning 
shell (taken at the layer of maximum energy production 
and at the onset of each pulse). Each mark locates a pulse. 
b Same as a, but for the maximum temperature reached 
at the base of the pulses. 



cleosynthesis calculations, presented in Appendix A, show 
that the abundance of 14 N emerging from H-burning ac- 
tually depends slightly on the temperature at which hy- 
drogen burns. Taking those temperature-dependent 14 N 
yields, the expected X£>( 22 Ne)/Xi„ i4 ( 22 Ne) ratio ranges 
between 149 and 151 (dotted line in Fig. ||a). Compar- 
ing that value with the intershell 22 Ne abundance pre- 
dicted by the stellar models (filled and open symbols in 
Fig. ||), we see that essentially all 14 N transforms into 
22 Ne at pulse temperatures T p (given in Fig. ||b) exceed- 
ing ~ 210 x 10 6 K. Such temperatures are already reached 
in the very first pulses of most of our model stars. The 
assumption of a complete conversion of 14 N into 22 Ne in 
the intershell layers is thus very good. 

At T p exceeding 250 - 280 x 10 6 K, on the other hand, 
Fig. H shows that the 22 Ne abundances predicted by the 
standard models are lower than that expected from the 
total conversion of the 14 N left over by the HBS. This is 
explained by the partial destruction of 22 Ne, which begins 
to burn by a capture. The 22 Ne destruction factor, how- 
ever, is less than 5% in all our stars, except in the 4 and 
6 M Q stars where it respectively reaches 10 and 20% in the 
24th pulse (with T p = 312 and 320 x 10 6 K, respectively). 
Such a destruction does not have significant implication 
on the 23 Na predictions for the 4M Q star, though, and is 
neglected in the synthetic and analytical calculations pre- 
sented in Sect. j|. The case of the 6 M Q star, on the other 
hand, is discussed in Sect. |5|. 



4. Sodium abundance predictions in models with 
dredge-up 

In this section, the production of sodium is analyzed in 
those AGB stars which do not experience HBB. The case 
of stars experiencing HBB is considered in Sect. [5| 



4-1- Synthetic calculations: the method 

In order to predict the surface abundance evolution of 
12 C, 22 Ne and 23 Na in our model stars, we have to use 
'synthetic' calculations since standard models cann ot re - 
produce consistently the 3DUP phenomenon (Sect. |3.l[) . 
In such synthetic calculations, the chemical evolution of 
a star is computed by estimating a posteriori the effects 
of successive 3DUP episodes on the intershell and surface 
abundances. In doing so, we assume that the overall struc- 
tural evolution of the model stars is not affected by the 
operation of those 3DUPs. The evolution, as a function of 
pulse number, of the mass M env of the envelope (Fig. 10 of 
MM99), of the masses M A , M B and M D of zones A, B and 
D, respectively (Fig. 8 of MM99), and of the mass fraction 
12 X l D of 12 C emerging from each pulse (Fig. [|), are thus 
taken from the standard stellar model predictions. The 
chemical evolution of the intershell layers and of the en- 
velope are then computed by considering their abundance 
changes during each interpulse/pulse/afterpulse phases. 

Let us consider pulse i. The abundances of 14 N and 
23 Na available in zone B ( 14 X B and 23 X B , respectively) 
are determined by the nucleosynthesis occurring in the 
HBS during the preceding interpulse i — 1. We assume that 
all available 12 C and 16 are converted in the HBS into 
14 N through the CNO cycles, and that all available 22 Ne 
is burned into 23 Na by the NeNa chain. The destruction 
of 23 Na by p-captures (at T > 60 x 10 6 K) is taken into 
account by a destruction factor d = A( 23 Na)/[A ( 23 Na) + 
||Xo( 22 Ne)] as defined in Appendix A. Its value is deter- 
mined in that Appendix as a function of temperature at 
which H burns, taken here to be the temperature of the 
HBS predicted by the standard model predictions (Fig. |5J). 
The slight contribution from 20 Ne burning, on the other 
hand, is neglected since 20 Ne is not produced in the inter- 
shell layers. This assumption imposes a maximum value 
of d — 1 . We then have 



14 



14 



xh = -; x 12 a 



14 



12 



2:S 



23 



X B =d[^x 22 X\ 



xir 



2:S 



— x w X l 
16 X Xs 



XI 



(1) 



(2) 



where 12 A|~ 1 , 22 AJ~ 4 and 23 A] _1 are respectively the 
abundances of 12 C, 22 Ne and 23 Na in the envelope dur- 
ing interpulse i — 1. 

During pulse i, we assume that all 14 N available in zone 
B is transformed into 22 Ne. In doing so, we neglect the fact 
that the transformation is incomplete in the coolest first 
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Table 3. Surface mass fractions of the 12 C, 14 N, 16 0, 22 Ne and 23 Na predicted by the standard model stars at the 
onset of the TP-AGB phase. The total C+N+O mass fraction is also given in the 4th line. The increase in those mass 
fractions relative to their initial values, due to the operation of the first and, possibly, second dredge-ups, are given in 
parentheses (in dex). 





4M 
Z=0.018 


3M 
Z=0.020 


2.5 M Q 
Z=0.018 


1.5 M Q 
Z=0.018 


2.5 M 
Z=0.008 


2.5 M 
Z=0.004 


12p 

14 N 

16 

CNO 

22 Ne 
23 Na 


1.785e-3 (-0.21) 
2.661e-3 (+0.40) 
8.632e-3 (-0.02) 
1.308e-2 (+0.00) 

1.071e-4 (-0.06) 
4.883e-5 (+0.19) 


2.013e-3 (-0.20) 
2.823e-3 (+0.38) 
9.699e-3 (-0.02) 
1.454e-2 (+0.00) 

1.212e-4 (-0.05) 
5.200e-5 (+0.17) 


1.820e-3 (-0.20) 
2.421e-3 (+0.36) 
8.849e-3 (-0.01) 
1.309e-2 (+0.00) 

1.105e-4 (-0.05) 
4.536e-5 (+0.16) 


2.052e-3 (-0.15) 
1.947e-3 (+0.27) 
9.100e-3 (-0.00) 
1.310e-2 (+0.00) 

1.216e-4 (-0.01) 
3.370e-5 (+0.03) 


7.714e-4 (-0.22) 
1.166e-3 (+0.40) 
3.882e-3 (-0.02) 
5.819e-3 (+0.00) 

4.760e-5 (-0.06) 
2.171e-5 (+0.19) 


3.660e-4 (-0.24) 
6.181e-4 (+0.42) 
1.926e-3 (-0.02) 
2.910e-3 (+0.00) 

2.311e-5 (-0.08) 
1.156e-5 (+0.22) 



pulses of some model stars. We further assume that the 
amount of 22 Ne burned in the hottest pulses of our stars 



is negligible (see Sect. 3.4). Knowing that sodium remains 
unaffected by He-burning in the pulse, the abundances of 
22 Ne and 23 Na emerging in zone D are respectively given 

by 



22^i _ 22 14 j M B 

A n — T - : X A f 



14 



M 



D 



Ji D — 



23 X B x M B + 23 Xi7 1 x (M l D - M B ) 



Ml 



(3) 



(4) 



The 3DUP process during afterpulse i then mixes part 
of the material of zone D into the envelope. The mass 
Mdup of material dredged-up from zone D is usually ex- 
pressed relative to the core mass increase AM C during 
the preceding interpulse by the 'dredge-up efficiency' pa- 
rameter A = Mdup I 'AM c . A value of A = 1, for example, 
means that the amount of mass dredged-up is equal to 
the amount of H-rich material burned during the inter- 
pulse (the core mass would then remain constant from 
one pulse to the next). Since AM C — Mb, the surface 
abundance of, for example, 23 Na after dredge-up is given 

by 



6 x l = 



23 Xr 1 x Mj nv + 23 Xh x\xM B 
ML,, + A x Mi 



(5) 



Similar relations hold for the surface abundances of 12 C 
and 22 Ne. 

Dredge-up is assumed to operate from the 5th pulse 
on in all the stars, and Eqs. ^ to (||) are solved pulse 
after pulse for the abundances in zones B, D and in the 
envelope. M env , Mb, Md and 12 Xo are taken from the 
standard AGB model predictions. We still need to specify 
the surface abundances at the onset of the thermally puls- 
ing AGB (TP-AGB) phase and the dredge-up efficiency. 
The abundances at the onset of the TP-AGB phase are 
taken from the standard model predictions as summarized 
in Table y. The dredge- up efficiency, on the other hand, is 
a function of stellar mass, metallicity and pulse number. 
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Fig. 6. Surface 23 Na versus 12 C mass fraction (in loga- 
rithm) predicted by synthetic calculations (Sect. |4.2| ), for 
the models as labeled in the figure. The abundances ex- 
pressed in dex are also shown in the upper and right axis, 
with [X]=log(X/X ). The 4M , Z = 0.018 synthetic star 
shows little variation in its surface abundances (at most 
0.3 dex variation in [ 12 C] between the first and 24th pulse 
with the parameters used in the synthetic calculations), 
and is not shown in the figure for clarity. 



Its value however is still a matter of research. Fortunately, 
the exact value of A is not crucial for our purpose, as long 
as the surface 23 Na abundance is expressed as a function 



of the surface 12 C (or C+N+O) abundance (see Sect. 4.3) 
For simplicity, we thus take A = 1. 



4-2. Sodium abundance predictions 

The surface 23 Na abundance evolution predicted by the 

are shown in 



synthetic calculations outlined in Sect. 4.1 
Fig. las a function of the surface 12 C abundance (each 
mark in the figure locates a pulse) . The carbon abundance 
is seen to increase steadily with pulse number, as expected 
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from the operation of the 3DUP episodes. The sodium 
abundance, on the other hand, remains unaffected during 
the first pulses, but steadily increases at more advanced 
pulses. 

We can roughly estimate the increase in the 23 Na 
surface abundance (A 23 A,,.) expected to result from 
given CNO abundances in the envelope. As explained in 



Sect. 2.2, the process of surface 23 Na enhancement from 



the envelope's CNO abundances requires two consecutive 
3DUPs, one to dredge-up 22 Ne synthesized in the H+HeBS 
from the CNO elements, and one to dredge-up 23 Na syn- 
thesized from 22 Ne during the following interpulse. At first 
guess, the resulting increase in the surface sodium abun- 
dance would thus be 



A 23 A 



j a 23 
d x f x — x 

14 



CNO 



Y„. 



(6) 



where / = Mg up /M env is the dilution factor suffered by 
the C-rich material dred ged- up in the envelope, d the 23 Na 
destruction factor (Sect. 4.1), and 



cno v _ if 12xr 

r s _ 12 a s 



4 A S + ^ 16 A S 
16 



(7) 



is the equivalent 14 N mass fraction of the CNO abundances 
in the envelope (i.e. the 14 N mass fraction which would re- 
sult from the total transformation of the CNO elements 
into 14 N). 12 A S , X s and 1G X S are respectively the surface 
mass fractions of 12 C, 14 N and 16 0. Taking, for example, 
d = 1, / = 0.005 and solar CNO abundances, one gets 
A 23 A S ~ 5 x 10~ 7 . This is 60 times lower than the solar 
sodium abundance, so that no significant enhancement in 
the surface 23 Na abundance would be predicted. This con- 
clusion holds true at all metallicities if we assume an initial 
solar distribution of the elements heavier than helium. 

The process of 23 Na production in the envelope is, how- 
ever, non linear, and Eq. (|6J) becomes invalid after the 
operation of few 3DUPs. This results from the fact that 
22 Ne builds up in the envelope dredge-up after dredge-up. 
Since 22 Ne is the seed nucleus for the production of 23 Na, 
the surface 23 Na abundance 23 A S would depend quadrati- 
cally on the surface 22 Ne abundance 22 X S (a constant 22 X S 
would lead to a linear increase of 23 X S with pulse num- 
ber). If we further take into account the increase in the 
envelope's carbon mass fraction pulse after pulse due to 
the dredge-up of primary 12 C, then 23 X S becomes roughly 
proportional to ( 12 X S ) 3 (see Appendix B). This explains 
the rapid increase of the surface 23 Na abundance as the 
surface carbon abundance increases. 

Figure O further reveals that the surface sodium abun- 
dance enhancements are higher in stars with lower metal- 
licities. This is due to the primary nature of 23 Na dredged- 
up to the surface. Actually, the [ 12 C] — [ 23 Na] predictions 
become independent of metallicity after the operation of a 
sufficient number of 3DUPs. For example, the surface mass 
fraction of 23 Na reaches ~ 10 -4 after about 25 3DUPs in 
the 2.5 M Q stars regardless of their metallicity (Fig. g). 
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Fig. 7. Same as Fig.g, but for a 2.5 M© star at differ- 
ent metallicities. The predictions at Z = 0.018 , 0. 008 and 
0.004 are from synthetic calculations (Sect. |4.2| ), while 
those at Z — 0.0 01 and 0.0001 are from analytical evalu- 
ations (Sect. 4.3). 



4-3. Discussion 

It is possible to further study the production of 23 Na in 
AGB stars analytically if we assume the constancy of the 
following four parameters: the intcrshell 12 C mass frac- 
tion 12 Xd dredged-up to the surface during each 3DUP 
episode, the dilution fraction / suffered by the C-rich ma- 
terial when dredged-up from the intershell layers into the 
envelope, the dilution factor h = Mb/Md suffered by the 
ashes of HBS ( 14 N in particular) when mixed from zone B 
of the intershell layers into the pulse, and the destruction 
factor d of 23 Na by p-captures. The analytical derivation is 
performed in Appendix B. It enables, in particular, to dis- 
entangle the effects of each of the fou r pa rameters, to ana- 
lyze t he in fluence of metallicity (Sect. |4.4| ) and stellar mass 
(Sect. 4.5 ) on the sodium predictions, and to e stim ate the 
uncertainties affecting those predictions (Sect. |4.6| ). 

The two main conclusions of the analytical study per- 
formed in Appendix B are: 

a) primary 23 Na production in zero metallicity starsn fol- 
lows a specific track in the [C+N+O] — [ 23 Na] diagram 
(with [C+N+0] = [ CAro r s ]), given by Eq. @ in the Ap- 
pendix and represented by a long-dashed line in Fig. H. It 

5 The evolution of Z = stars might qualitatively be quite 
different than that of Z 7^ stars. However, the only require- 
ment for our sodium production scenario to work is the occur- 
rence of 3DUPs which bring primary carbon into the envelope. 
And model predictions confirm the existence of thermal insta- 
bilities in the He-burning shell of AGB stars - and thus most 
probably of 3DUP events - at me tallicities down to Z=0 for 
M ^ 4M (Fujimoto et al. |l984p . In any case, the LOPSE 
could very well be defined by the asymptotic evolution of 23 X S 
in a very low (non zero) metallicity star. 
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is called the line of primary sodium enrichment (LOPSE), 
and is, to first order, a straight line of slope 3 [cf. Eq. (pi) 
in the Appendix]; 

b) in non-zero metallicity stars, the surface 23 Na abun- 
dance is given by Eq. (16) in the Appendix. 23 X S is pre- 



dicted to keep constant during the first 3DUPs until the 
abundance of primary 23 Na produced by the star (and 
given by the LOPSE) reaches a level comparable to the 
initial surface 23 Na abundance. It then increases steadily, 
and approaches the LOPSE asymptotically. 

The evolutionary tracks predicted at Z = 10~ 4 and 
1(T 3 in the [C+N+0]-[ 23 Na] diagram by Eq. © of the 
Appendix, for example, are shown by the two unmarked 
solid lines in Fig. 0. Those analytical predictions are in 
good agreement with the 23 X S evolution predicted by the 
synthetic calculations. For comparison, the synthetic pre- 
dictions for the 2.5 M Q stars at Z = 0.004, 0.008 and 0.018 
are also drawn in Fig. (symbol-marked solid lines); 

Appendix B further analyses the sensitivity of the an- 
alytical 23 Na abundance predictions to 12 Aq, /, h and d. 
The value of 12 Xrj predicted by our standard models has 



been analyzed in Sect. 3.4. It reaches a maximum value 



of 25-30% (depending on the initial stellar mass) after a 
dozen of pulses, and asymptotically decreases thereafter 
to a 'canonical' value of 12 X D = 0.20 (Fig. ||). This value 
is the one adopted in the curves at Z = 10~ 4 and 10 -3 
displayed in Fig. |7| The analysis of Appendix B reveals a 
strong sensitivity of the LOPSE on 12 Xd, though. This re- 
sults from a dependence of 23 X S on the square of 12 Xrj [cf. 
Eq. ( |l9| ) in the Appendix] . Increasing 12 Xrj by a factor of 
two, from 0.2 to 0.4 for example, decreases 23 X S predicted 
by the LOPSE at a given surface C+N+O abundance by 
a factor of about four (0.6 dex; Fig. [lja of Appendix B). 
Fortunately, the value of 12 Xu after a dozen of pulses is 
rather well constraint according to standard stellar mod- 
els, lying between 0.20 and 0.24. Furthermore, the valu e 
of 12 Xr> is rather independent of metallicity (Sect. 3.4). 



Assuming a constant value for that parameter is thus a 
good enough assumption for our purposes. Future evolu- 
tionary AGB model calculations consistently reproducing 
the 3DUP processes should lead to more refined predic- 
tions. 

The second parameter is the dredge-up dilution fac- 
tor /. Appendix B reveals that 23 X S predictions are not 
very sensitive to /. To first order, the LOPSE is even in- 
dependent on that parameter [Eq. ( |l9| ) in the Appendix] . 
This independence of 12 X S on / is welcomed since the 
dredge-up efficiency as a function of mass and metallic- 
ity is still a matter of research. We choose a 'canonical' 
value / = 0.005. This value corresponds to that predicted 
by the standard 2.5 Mq models if we assume a dredge-up 
efficiency of A = 1, as shown in Fig. ||a. 

The third parameter is dilution factor h of H-burning 
ashes into pulses. Standard AGB models predict a contin- 
uous increase of that parameter along the AGB, from 0.05 
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Fig. 8. a Dilution factor / = Md U p/M env , as a function 
of core mass, suffered by the material dredged-up by a 
3DUP episode from the C-rich intershell layers into the 
envelope, as predicted by the standard AGB models when 
a dredge-up efficiency of A = 1 is assumed. The symbols 
identifying the mass and metallicity of each line are the 
same as those Fig. [|. b Same as a, but for the dilution 
factor h = Mb /Mo suffered by the ashes of the hydrogen 
burning shell when mixed into the pulses. 



to 0.50-0.60 as reported in Fig. pp. We use a 'canonical' 
value of h = 0.40. The analysis developed in Appendix B, 
however, reveals a linear dependence of 23 X S to h on the 
LOPSE [Eq. (19) in the Appendix]. As a result, an in- 
crease of h from 0.2 to 0.4, for example, would result in 
an increase of 23 X S by a factor of two at a given 12 X S . Giv- 
ing a realistic range 0.20-0.60 for h in AGB models, the 
error on 23 X S due to h thus amounts to a factor of three. 
The last parameter is the 23 Na destruction factor d. 
With the nuclear reaction rates used in the standard 
model calculations, d drops below 1 only for Tjj > 70 x 
10 6 K (cf. Appendix A). Because Tjj is predicted to be 
lower than that value in most of our standard model stars 
(Fig. ||), we adopt a 'canonical' value of 1. Since 23 X S de- 
pends linearly on d, the error introduced by using d = 1 



is of a factor of three (see, however, Sect. 4.6). 

The analytical equation (Ufl) of the Appendix thus of- 
fers a powerful tool to analyze the production of primary 
23 Na in AGB stars. The main qualitative difference be- 
tween [C+N+O] — [ 23 Na] predictions from the analytically- 
derived relation (uM in Appendix B and tho se from the 
synthetic calculations presented in Sect. L2 (and which 
take into account the full stellar characteristics as pre- 
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dieted by standard AGB models pulse after pulse) is the 
initial 23 Na abundance at the onset of the TP-AGB phase 
(see Fig. R). The analytical calculations begin with an ini- 
tial solar "Na abundance scaled to the metallicity (dotted 
line in Fig. yy, while the synthetic calculations take into 
account the Na abundance enhancements resulting from 
first and second dredge-ups. The increase in 23 X S due to 
first dredge-up is, however, less than 0.3 dex in the stars of 
interest (cf. Table pf). Its contribution becomes meaning- 
less when the star reaches the LOPSE, i.e. when primary 
23 Na dominates sodium production. 

4-4- Sodium production as a function of metallicity 

Figure ffl shows that the increase in the surface C+N+O 
abundance required to reach the LOPSE is greater in lower 
metallicity stars. At solar metallicity, for example, the 
surface C+N+O abundance must increase by less than 
0.4 dex before it reaches the LOPSE, while 1.2 dex are re- 
quired at Z = 10~ 3 , and 1.9 dex at Z = 10~ 4 . Yet, metal- 
poor stars are expected to be more efficient in produc- 
ing primary 23 Na than solar-metallicity stars. This results 
from the fact that the initial surface C+N+O abundance 
is proportional to the metallicity (the CNO elements rep- 
resent 72% of the metallicity), while the primary 12 C mass 
fraction produced in the HeBS is independent of metallic- 
ity (Fig. |). As a result, the relative increase of the surface 
12 C abundance is higher in metal-poor stars. We can easily 
estimate the number n^up of 3DUPs necessary to increase 
the surface C+N+O abundance by a factor of, let us say, 
x. Let us assume for simplicity that the dredge-up dilution 
factor / is constant from one pulse to the next. We then 
have 



n du p = (*-!) 



0.72 x Z 
12 X D xf 



(8) 



For example, the number of 3DUPs necessary to increase 
'' NO Y s by a factor of three (i.e. by 0.5 dex) is, assuming 



12 



X D = 0.24 and / = 0.006, 



Tldup 



1000 z. 



This corresponds to 18 dredge-ups in a solar metallic- 
ity star, 8 at Z = 0.008, and only one at Z < 0.001. 
The efficiency of surface C enrichment thus increases very 
rapidly with decreasing metallicity. Low-metallicity stars 
should thus be good candidates for the production of pri- 
mary 23 Na. 

At Z — 10~ 4 , for example, the number of 3DUPs re- 
quired for CNO Y s to reach the LOPSE (i.e. to increase by 
1.9 dex) is expected to be 4 according to Eq. (||). Of course, 
Eq. (g) only provides an order of estimate, and those pre- 
dictions should be confirmed by future evolutionary AGB 
models consistently reproducing the 3DUPs. In particular, 
the dependence law of the dredge-up efficiency on mass, 
metallicity and pulse number is of crucial importance in 
determining the number of pulses required to reach the 



LOPSE. According to current estimations the dredge-up 
efficiency is expected to increase with decreasing metallic- 
ity. This strengthens our conclusion that low-metallicity 
stars should be efficient sites for the production of pri- 
mary 23 Na. 

4-5. Sodium production as a function of stellar mass 



The synthetic calculations presented in Sect. 4.2 reveal a 



dependence of the sodium abundances on the stellar mass. 
Figure ^ indeed shows that 23 X S begins to increase at a 
lower C+N+O abundance in more massive stars. In other 
words, 23 Na production would be favored in more massive 
AGB stars. 

Why such a trend in the synthetic predictions? The 
analytical study performed in Appendix B again offers 
the necessary tool to answer that question. What actually 
happens is that in a massive AGB star, the low dilution 
factor / (M env being high) leads to a slow increase of 
the surface 12 C abundance with pulse number. It takes 
then many 3DUP episodes for the surface C+N+O abun- 
dance to reach the LOPSE. By that time, the h factor 
has already increased to a high value according to stan- 
dard model predictions (Fig. ^b), leading to higher 23 Na 
abundances. In a low-mass AGB star, on the contrary, the 
high dilution factors / lead to a rapid increase in the sur- 
face 12 C abundance with pulse number. As a result, h is 
still low when the surface C+N+O abundance reaches the 
LOPSE, leading to lower 23 X S values at a given C+N+O 
surface abundance. 

The 23 A S predictions thus also depend on the time at 
which 3DUP becomes efficient on the TP-AGB (thereby 
determining h), i.e. on the history of the 3DUP efficiency. 
In the synthetic calculations presented in Sect. 



4.1 



it is 

assumed that 3DUP begins to operate at the 5th pulse 
irrespective of the stellar mass and metallicity, and that 
its efficiency is equal to A = 1 for all subsequent pulses. In 
reality, the dredge-up efficiency is expected to vary with 
stellar mass, metallicity, pulse number and mass loss rate. 
The lack of reliable predictions, at the present time, for 
the 3DUP efficiency as a function of pulse number thus in- 
troduces an additional source of uncertainty on the quan- 
titative prediction of 23 A S . 

4-6. Uncertainties on the 23 Na abundance predictions 



The quantitative prediction of primary 23 Na prod ucti on 
in AGB stars has been shown in Sects. 4.3 and 4.5 to 
depend sensitively on the 12 C abundance produced in the 
HeBS, on the intershell structural parameter h, and on 
the nuclear destruction factor d. 

In our standard AGB models, the intershell 12 C abun- 
dance reaches the asymptotic value of 0.20-0.24 after a suf- 



ficient number of pulses (Sect. 3.4). 'Non standard' input 



physics, ho wever , may lead to a different conclusion. Her- 
wig et al. ( 1997 ), for example, obtain l2 Xry = 0.50 when 
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applying extra- mixing (overshooting) at the borders of the 
convective pulses. Since the amount of overshooting at the 
convective pulse borders is certainly not yet established, 
we have to consider 12 J*Cd to be known to only a factor of 
2.5. Moreover, an additional source of discomfort in the 
quantitative prediction of 12 Ad results from the uncer- 
tainty affecting the nuclear reaction rate of 12 C (a , 7) 16 



(see Arnould et al. 1999). The role of Ad is all the more 



important since the primary 23 Na abundance predicted by 
the LOPSE depends on the square of 12 Xo- If we thus take 
a global uncertainty factor of about three on 12 Xp due to 
convection and nuclear reaction rates, then the sodium 
abundance predicted at a given C+N+O surface abun- 
dance is reliable only to a factor of about 10. 

The uncertainties on h, on the other hand, have been 



shown in Sect. 4.3 to lead to an uncertainty factor on 23 X S 
of about 3. 

Finally, the uncertainty on d has been shown in 



Sect. 4.3 to lead to an uncertainty on 23 X S of a fac- 



tor of 3. The new NACRE compilation, however, predict 
much higher 23 Na destruction rates than before (Sect. |373J). 
But the very large uncertainties still affecting the p- 
capture rates on 23 Na prevent from giving a definite con- 
clusion about their impact on 23 A S predictions. A 23 Na 
destruction factor of 'only' 3, for example, is still within 
the NACRE uncertainties. We thus adopt this optimistic 
value, keeping however in mind a possible higher destruc- 
tion factor due to nuclear reaction rates. 

All together, then, the present uncertainty affecting 
23 A S predictions at a given surface C+N+O abundance 
amounts to two orders of magnitude. 

4-7. Summary 

Both the predictions of synthetic AGB calculations pre- 



sented in Sect. 4.2 and the analytic considerations of 



Sect. 4.3 do support the scenario described in Sect. 2.2 for 



the production of primary sodium in AGB stars. The evo- 
lution of the surface 23 Na abundance obeys a well defined 
pattern as a function of the surface C+N+O abundance, 
given by Eq. ( |lq ) in Appendix B and shown in Fig. M for 
several metallicities. 

All AGB stars are thus potential sites for the produc- 
tion of primary 23 Na. For such a production to be effi- 
cient, a sufficient number of 3DUPs must occur in order 
to increase the surface C+N+O abundance to the level 
required by the LOPSE (dashed line in Fig. pp at the 
given initial 23 A S . The amount of primary sodium pro- 
duced in a given AGB star thus depends on the num- 
ber of 3DUP episodes experienced by that star. Since 
the dredge-up efficiency increases with decreasing metal- 
licity, low-metallicity stars are expected to be more effi- 
cient in producing primary 23 Na. This conclusion is fur- 
thermore strengthened by the fact that the relative in- 
crease in the surface carbon abundance during a 3DUP 
episode is higher in low than in high metallicity stars. 
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Fig. 9. a Evolution of the temperature at the base of the 
envelope of a 6 M©, Z = 0.018 standard model star during 
the first 30 pulses, b Same as a, but for the surface 23 Na 
mass fraction. The origin of time is set to the first pulse. 

The uncertainties affecting both synthetic and analytic 
predictions of 23 X S are, however, still very large at present 
time (two orders of magnitude). In particular, the uncer- 
tainties affecting the 23 Na burning reaction rates is of great 
concern. Eventual observation of 23 Na overabundance in 
AGB stars and related objects would provide a definite 
proof of the feasibility of our scenario for the production of 
primary 23 Na. This question is briefly addressed in Sect, pi 
Before doing that, however, the case of massive AGB stars 
experiencing HBB is discussed in the next section. 

5. Sodium production by hot bottom burning 

In massive AGB stars (M £ 5 M Q at Z=0.02), the temper- 
ature at the base of the envelope may reach high enough 
values for H-burning to operate in the envelope. This 
process, called HBB, alters the surface abundances of H- 
burning nuclei without the need to invoke any dredge-up 
process. In particular, the activation of the NeNa chain 
of H-burning naturally leads to the production of 23 Na 
when the temperature at the base of the envelope exceeds 
60 x 10 6 K. 

5.1. Secondary sodium 

In the absence of any 3DUP event, secondary sodium can 
be synthesized by HBB from the initial 22 Ne present in the 
envelope. The evolution of the temperature at the base of 
the envelope of the 6M©, Z = 0.018 standard model star, 
for example, and of its surface 23 Na abundance are shown 
in Fig. |§| during the first 30 pulses. The production of 
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secondary sodium is seen to become efficient when the base 
temperature increases above 70 — 80 x 10 6 K. The increase 
in the surface 23 Na abundance is steady, eventhough not 
very rapid (30% increase over about 25 pulses). 

5.2. Primary sodium 

In the presence of 3DUPs, the increase in the surface 23 Na 
abundance through HBB may be considerably higher due 
to the dredge-up of primary 22 Ne. In such a scenario, 
only one dredge-up/interpulse/pulse/dredge-up event is 
required to dredge-up 12 C and process it to 22 Ne, the con- 
version of 22 Ne to 23 Na being subsequently performed in 
the envelope. The fact that the dredge-up efficiency in- 
creases with stellar mass renders this scenario very at- 
tractive. Several important points must however be con- 
sidered. 

First, the standard models reveal a destruction of the 
intershell 22 Ne by a-capture at the He-burning shell tem- 
peratures characterizing massive AGB stars (Sect. 3.4). 
In the 6M Q model star, the destruction factor amounts 
to 20% at the 24th pulse (Fig. ph. Let us, however, neglect 
this destruction factor for now. 

Second, the dilution factor / of the 22 Ne-rich mate- 
rial mixed from the intershell layers into the envelope is 
predicted by standard models to be only 0.0002, much 
smaller than in lower mass stars (/ ~ 0.001, 0.003, 0.005 
and 0.01 in the 4, 3, 2.5 and 1.5 M© model stars, respec- 
tively, see Fig. @). As a result, the surface 22 Ne mass frac- 
tion is not expected to increase very rapidly as a result of 
3DUP episodes in massive AGB stars. For example, a syn- 
thetic calculation of the surface 22 Ne abundance evolution 
in the 6M© model star, performed with the recipe out- 
predicts only 20% increase in the surface 



4.1 



lined in Sect. 

22 Ne abundance after 24 pulses. The synthesis of primary 
sodium by 3DUP + HBB may thus not be very efficient. 
A more positive conclusion may however be obtained in 
low-metallicity stars, where the low initial 22 Ne abundance 
renders the relative contribution of the intershell 22 Nc to 
the envelope's 22 Ne abundance more important. 

Third, synthesizing primary sodium in the intershell 
layers of massive AGB stars by two successive dredge- 
up/interpulse/pulse sequences, in a similar way as in lower 
mass stars, is not efficient cither. The temperatures in the 
H-burning shell of those massive AGB stars indeed exceed 
75 x 10 6 K according to the predictions from the 6M Q 
model star (cf. Fig. ||) , and lead to the destruction of 23 Na 
by p-capture in the HBS. At the end of the 23rd interpulse 
of the 6 M Q model star, for example, the 23 Na intershell 
abundance is 30% lower than its surface abundance. 

Finally, we have to consider the effects of mass loss 
which characterizes AGB stars. It reduces the mass of the 
envelope, thereby decreasing the temperature at its base, 
and eventually switching HBB off. After how many inter- 
pulses is HBB switched off, however, depends on the still 
uncertain mass loss rates suffered by those stars. 



The production of primary sodium by HBB seems thus 
not favored in massive AGB stars. There are however 
still too many uncertainties, such as those concerning the 
dredge-up efficiency or mass loss rates, to allow a definite 
conclusion to be drawn. 

6. Observations 

The primary sodium, if produced efficiently in AGB stars, 
should reveal its imprint in the abundances at the surface 
of AGB and post- AGB stars. According to the analysis of 
Sect. [|, sodium could be produced in AGB stars of nearly 
solar metallicity when their 12 C abundance has increased 
by at least 0.7 dex (Fig. g). This would correspond to 
stars displaying C/O ratios above ~2. Unfortunately, the 
number of known such objects is low, due to the fact that 
those AGB stars are usually enshrouded into a circum- 
stellar dust which obscures the central AGB star. To my 
knowledge, no sodium abundance has so far been mea- 
sured in those objects. 

The analysis of post-AGB stars seems more promis- 
ing. Those objects are expected to display the abundance 
patterns characterizing AGB stars at the end of their life, 
and in particular high C+N+O abundances if 3DUP oc- 
curred during their AGB phase. Spectroscopic analysis 
of some of those objects has been performed during the 



last decade (e.g., Parthasarathy et al. 1992, Gonzalez & 



Wallcrstcin |1992j Zacs et al. |199§ Van Winckel et al. |1996 , 
Decin et al. 199S). In particular, Gonzalez & Wallcrstcin 



(1992) analyze the chemical composition of a post-AGB 



star, ROA 24, belonging to the globular cluster w Cen. 
They find significant overabundances of C, N, O, and s- 
process elements compared to the abundances measured in 
other giants in u> Cen. The fact that the s-process elements 
are enhanced is a very strong indication of the post-AGB 
nature of that star. Very interestingly, they also find a 
significant overabundance of Na in that object relative to 
the abundance in other giants of u Cen. This would be the 
first observational evidence of primary 23 Na production in 
AGB stars. A more thorough analysis of that case and of 
other post-AGB stars is outside the scope of this article, 
though, and is the object of a separate study (Mowlavi & 
Van Winckel, in preparation). 

7. Conclusions 

The main conclusion of this paper is that AGB stars are 
potential sites for the production of primary sodium. For 
the scenario presented in Sect. @ to be efficient, however, 
a certain level of surface C+N+O abundance must be 
reached by the operation of sufficient number of 3DUP 
events (Sect. §). 

The surface 23 Na abundance evolution can be pre- 
dicted as a function of the surface C+N+O abundance 
by the analytical relation ( |l6[) presented in Appendix B, 
with 'recommended' values of 12 Ad = 0.20, / = 1, h = 0.4 
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and d = 1. The values of 12 Xd, h and d are suggested from 
standard AGB model predictions, while the value of / is 
suggested after the study of Mowlavi ( |1999[ ). The 23 Na 
abundance predictions are sensitive to 12 Ad , h and d, but 
not much on / as long as X s ( 23 Na) is expressed as a func- 
tion of the C+N+O abundance. 

An important result of the analytical study is the pre- 
diction of a LOPSE in the [C+N+O] -[ 23 Na] diagram, 
which provides the 23 Na abundance evolution in zero 
metallicity stars as a function of their surface C+N+O 
abundance. The surface 23 Na abundance in non-zero 
metallicity stars, on the other hand, is shown to evolve 
asymptotically towards the LOPSE. 

Hot bottom burning in massive AGB stars may fur- 
ther help the production of 23 Na by synthesizing it di- 
rectly in the envelope. In that case, only one dredge- 
up/interpulse/pulse/dredge-up sequence is necessary to 
synthesize and dredge-up primary 22 Ne, followed by its 
burning into 23 Na in the envelope. The analysis presented 
in Sect. |g, however, suggests that such primary sodium 
production may not be very efficient in those stars, at least 
at solar metallicity. Low metallicity stars may be better 
candidates for producing sodium through 3DUP + HBB. 

The uncertainty on the primary 23 Na abundance pre- 
diction at a given C+N+O abundance is presently of two 
orders of magnitude, even with synthetic AGB calcula- 
tions. This poor prediction partly results from our still 
poor quantitative knowledge of the 3DUP efficiency as 
a function of stellar mass, metallicity and pulse number, 
since the values of 12 Xd and h vary with pulse number in a 
given star, and partly from the uncertainties affecting the 
p-capture rates on 23 Na. Future AGB model predictions 
consistently reproducing the 3DUPs, and combined with 
improved reaction rates, would hopefully enable to provide 
a better quantitative prediction of sodium production in 
low- and intermediate-mass stars. 

Meanwhile, we can search for an observational evi- 
dence for this sodium production scenario in carbon stars 
and related objects. Such an evidence may already have 
been reported in at least one post- AGB star in the glob- 
ular cluster lu Cen (cf. Sect. ||). Further determinations 
of sodium abundance in those objects and, if possible, in 
planetary nebulae are eagerly called for. 

Finally, let me note that the efficient production of 
primary sodium in AGB stars, if confirmed by observa- 
tions, would imply a non negligible role of those objects 
in the chemical evolution of galaxies in general, and of our 
Galaxy in particular. 



Appendix A: The CNO and Ne-Na modes of hy- 
drogen burning 

The CNO cycles 
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Fig. 10. Mass fractions of the stable nuclides involved in 
the CNO cycles at H exhaustion [A(H) = IO -5 ] versus the 
temperature at which hydrogen burns. The dotted lines 
represent the yields using the NACRE reaction rates (see 
text). The initial (solar) abundances are indicated by the 
horizontal thick lines at the right of the figure. 



al. 1995 ). We only summarize here their main features in 
relation with the 23 Na production in AGB stars. The rates 
of the CNO reactions adopted in the st andar d model cal- 
culations are the same as in Mowlavi (1995), except for 



17 (p , a) 14 N and 17 (p , 7) 14 N whose rates are taken 



from Blackmon et al. ( 1995| ). 

The CNO yields obtained at hydrogen exhaustion in 
one-zone nucleosynthesis calculations with the tempera- 
ture assumed constant during the burning are displayed 
in Fig. 10 as a function of temperature. As is well known, 
essentially all CNO elements are transformed into 14 N. 
However, 12 C, 13 C and 16 are still present at H exhaus- 
tion at the level of few percents of the 14 N mass fraction. 
At temperatures > 40 x IO 6 K characterizing H-burning 
in AGB stars, the mass fraction of 12 C (which is the most 
abundant of the CNO elements after 14 N, see Fig. |10|) in- 
creases with temperature. Concomitantly, the 14 N mass 
fraction at H exhaustion decreases slightly with tempera- 
ture (few percents fr om 40 to over 70 x IO 6 K). This prop- 
erty is used in Sect. 3.4 to analyze the 22 Ne production in 
AGB stars. 

The impact of the new NACRE reaction rates 
(Arnould et al. 1999) on the CNO yields of hydrogen burn- 
ing is shown by the dotted lines in Fig. QQ. It is seen to be 
negligible for the N production. 

The NeNa chain 



The CNO modes of hydrogen burning are well known The reactions involved in the NeNa chain of hydrogen 
and documented in the literature (see, e.g., Arnould et burning are shown in Fig. [11]. The rates of p captures 
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Fig. 11. Reactions of the NeNa chain. Stable nuclides are 
enclosed in thick circles. The possible leakages in and out 
of the chain are represented by the dashed lines. 



on 20 Ne, 21 Ne, 22 Ne and 23 Na are taken from El Eid & 



Champagne (1995). 

The NeNa yields at hydrogen exhaustion are displayed 
in Fig. [l^ a (solid lines). The main features concerning 
23 Na production are the following (see Arnould et al. 1995 
for more details): 

- 22 Ne is entirely transformed into 23 Na at tempera- 
tures as low as 15 x 10 6 K. Consequently, the ratio d = 
A( 23 Na)/[A ( 23 Na) + §§A ( 22 Ne)], where X designates 
the initial (here solar) abundance, is equal to one up to 

- 35 x 10 6 K (solid line in Fig. ||b). 

- an extra production of 23 Na is predicted at temperatures 
exceeding 35 x 10 6 K, reaching - 60% at T ~ 60 x 10 6 K 
(Fig. |l2|b). This extra production results from the slight 
burning of the abundant 20 Ne; 

- 23 Na begins to burn by proton capture at T > 60 x 10 6 K, 
with d = 0.5 at T = 80 x 10 6 K. 

It is interesting to further note that the El Eid & 
Champagne (1995) predict 23 Na (p , a) 20 Ne reaction rates 
always faster than those of 23 Na (p , 7) 24 Mg , which en- 
sures that the NeNa chain is a cycle. 

The impact of the new NACRE reaction rates on the 
NeNa yields is shown by the dotted lines in Fig. p~2. The 
23 Na destruction by p-capture is seen to be much more 
important with NACRE than with the El Eid & Cham- 
pagne reaction rates above 55 x 10 6 K leading to about 
three times less 23 Na abundances at temperatures above 
65 x 10 6 K. Yet, the 23 Na abundance at hydrogen exhaus- 
tion does not decrease below its initial abundance at all 
temperatures of relevance for our study (i.e. d > rip for 
T < 80 x 10 6 K, see Fig. |Jb). Arnould et al. ( |l999| ) fur- 
ther note that the cycling character of the NeNa chain is 
put into question with the new NACRE rates, but still 
depends on the large uncertainties affecting the p-capture 
rates on 23 Na. 

It must be stressed that the large uncertainties affect- 
ing those p-capture rates render the 23 Na yields from the 
NeNa chain highly uncertain. A doptin g the most favorable 
NACRE rates, Arnould et al. ( 1999|) show, for example, 
that the 23 Na mass fraction at H exhaustion can still be 




40 60 

T (xlO 8 K) 

Fig. 12. a Same as Fig. nfl, but for the stable nuclides 
involved in the NeNa chain, b d = A( 23 Na)/[A ( 23 Na) + 
||Ao( 22 Ne)], where Xq designates the initial abundance. 
The solar value d© is indicated by the horizontal thick line 
at the right of the figure. 



10 4 at 80 x 10 6 K, i.e. similar to the results obtained with 
the El Eid & Champagne rates. 

Appendix B: Analytical study of the surface 23 Na 
abundance evolution in AGB stars 

In order to analytically derive the evolution of the sur- 
face 23 Na abundance as a function of the surface C+N+O 
abundance, we make the following simplifying assump- 
tions: 

(a) the intershell 12 C mass fraction 12 Ad left over by the 
pulses is constant from one pulse to the next; 

(b) the dilution factor / = Md U p/M env suffered by the 
intershell nuclei when dredged- up into the envelope during 
a 3DUP episode is constant from one 3DUP to the next; 

(c) the dilution factor h = Mb I Mb suffered by the inter- 
shell 14 N when injected from zone B (Fig. pi) into the next 
pulse is constant from one pulse to the next; 

(d) all 12 C and 16 are transformed into 14 N in the 
CNO cycles, and all 22 Ne into 23 Na in the NeNa chain. 
The extra-production of 23 Na from 20 Ne is neglected, 
and the 23 Na destruction by p-capture is taken into ac- 
count through the parameter d = A( 23 Na)/[A ( 23 Na) + 
||A ( 22 Ne)], where A designates the initial abundance. 
The parameter is assumed constant from one pulse to the 
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next. A value of h = 1 means no 23 Na destruction, while 
a destruction of 23 Na by p-captures leads to ft, < 1 ; 

(e) the intershell 14 N is assumed to be entirely transformed 
into 22 Ne in the pulses. The 22 Ne destruction by a-capture 



4.1) 



in the hottest pulses is neglected (cf. Sect. 

The vali dity of assumptions (a)-(d) is discussed in 
Sect. L3, and that of (e) in Sect. ^J. The sensitivity of 
the 23 Na abundance predictions to 22 Xd, /, h and d, on 
the other hand, is analyzed at the end of this Appendix. 
Let us denote the surface mass fractions of 12 C, 14 N, 
16 0, 22 Ne and 23 Na by c, n, o, e and s, respectively (i.e. 
c = 12 A S , ... , s = 23 X S ). Let us further denote by y 
the total abundance _ Y s of the CNO elements in the 
envelope (cf. Sect. 12). The evolution of the surface 12 C, 
22 Ne and 23 Na mass fractions, which are dictated by the 
3DUP episodes, are respectively given by 



f 12 x L 
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where the indice i (i — 1) refers to the abundances during 
interpulse i (i — 1), and 22 Xds and 23 Xr>,i are respectively 
the intershell abundances of 22 Ne and 23 Na in zone D af- 
ter pulse i. 22 Xd,i can be evaluated by recalling that the 
intershell 22 Ne abundance results from the transformation 
of the CNO elements into 14 N in the HBS, followed by the 
mixing of 14 N into the next pulse with a dilution factor h, 
followed by its transformation in the pulse into 22 Ne by 
two a-capture reactions. We thus have 
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y| Oi. Equation ( |1 0| ) then becomes 

(12) 

We note that the evolution of j/j is dictated by that of 
the surface carbon abundance, since the surface nitrogen 
and oxygen abundances are not affected by the 3DUP 
episodesH. We thus have 



14 14 

yi = 7T (k + n + — o 
12 lb 



(13) 



with a given by Eq. (g), and n and o being respectively 
the surface mass fractions of nitrogen and oxygen at the 
onset of the TP-AGB phase. 



Actually oxygen is also produced in the pulses by 
12 C (a , 7) 16 . However, its mass fraction in the intershell lay- 
ers is predicted by our standard models to be about 20 times 
less than the intershell 12 C mass fraction. For this reason, we 
neglect in this analytical study its contribution to the surface 
CNO abundance evolution. 



The intershell abundance 23 Arj of 23 Na, on the other 
hand, results from the conversion of 22 Ne into 23 Na by p- 
capture operating in the HBS. Since 23 Na is not destroyed 
in the pulses, we haven 
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and Eq. (|Tl]) becomes 

Si = Si_i +/ — fy-l. 

Equations @, (|T2|), ( |l3| ) and (|14| ) constitute a system 
of recursive equations in the variables Cj (or yi), ej and 
Si. It can be solved rather simply as a function of pulse 
number i. Knowing that 



'- 1 1 

i=i 
and 

EE k = jU(*-i)(*-2), 
3=1 fe=i 

the solution to Eqs. @, ©-(0) is 
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where cq, yo, eo and so are the surface abundances at the 
onset of the TP-AGB phase. 

An analytical relation for the surface 23 Na abundance 
as a function of that of C+N+O can now easily be ob- 
tained by eliminating the variable i between Eqs. (15) 
for Ci and Sj. We have, however, to be cautious that two 
consecutive pulses are necessary before 23 Na can build up 
from the envelopes CNO elements. We thus take the sur- 
face abundances C2, yi, &i and S2 at the second pulse as 
initial conditions in Eqs. (15) instead of the surface abun- 
dances at the onset of the TP-AGB phase [replace Co, yo, 
eo, so respectively by C2, 2/2, &2, S2 in Eqs. (|l5[, and elim- 
inate i between the expressions for a and s, so obtained. 
This leads to a relation between s and y valid from the sec- 
ond 3DUP on, i.e. for Vi > y + 2 f 12 X D . The surface 23 Na 
abundance after the first 3DUP (i.e. at y = yo + f 12 Xjj), 
on the other hand, is immediately obtained from Eqs. (|15| 



In reality, of course, we have 23 Xn,i = dx[Si 



1 23 
-l+oo 



For simplicity in the analytical derivations, however, we ne- 
glect the contribution of Si-i compared to that of d-i. This 
approximation leads to a slight underestimation of s; . 
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by taking i = 1. The result is, after some elementary al- 
gebraic manipulations, 
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with [from Eqs. fll3j) and (|15|)] 
y 2 = yo + 2 / ^ 12 X D 
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Zero metallicity stars 

Let us first consider the 23 Na production in zero metal- 
licity stars]] (i.e. the production of pure primary sodium) . 
Putting Co = «o = °o = eo = So = in Eqs. ( ]l6|) and (|l7J), 
we get 



a = ^ d ft ^ 

12 



3 ' 



2/ 12 X D 
12 X D 

c-2f 12 X D 



\f 



-2f 12 X D Y 
12 X D ) 



12 



■Xd 



(18) 



Equation ([18]) is shown in thick line in Fig. [13] for 
'canonical' values of the parameters 12 Xp = 0.2, / = 
0.005, h = 0.4 and d — 1 (see Sect. 4.3). It represents 
the abundance of primary 23 Na resulting from the sole 
contribution of the primary 12 C produced in the HeBS. It 
is therefore called the line of primary sodium enrichment 
(LOPSE). 

The sensitivity of the LOPSE to the parameters 12 Xd, 
/, h and d can be estimated from Eq. (Uq). Knowing that 
/ << 1 in most circumstances, we can neglect the terms 
in / and f 2 in the right hand side of Eq. ( |l8| ) , and write 
to first order 
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Fig. 13. Surface 23 Na abundance, expressed in dex rela- 
tive to the solar abundance {i.e. log[X( 23 Na)/X ( 23 Na)]}, 
predicted analytically by Eq. ( J16[ ) for different metallici- 
ties Z as labeled in the figure. The predictions are given as 
a function of the surface CNO abundances, also expressed 
in dex (by number) relative to the solar abundance. The 
'canonical' values used for the parameters 12 Xrj, f and h 
entering (Eq. [lq ) are given in the insert at the upper-left 
corner of the figure. The filled circles on each track locate 
the pulses as predicted by (Eq. 0) . The dashed line repre- 
sents the abundances of 23 Na versus the CNO abundances 
if a solar distribution of elements heavier than helium is 
assumed at different metallicities. The thick line is the 'line 
of primary 23 Na production' defining the 23 Na abundance 
evolution in zero metallicity stars (see text). 



Several important conclusions can immediately be 
drawn from Eq. (]l9|): 

1 . the surface 23 Na abundance increases as the third power 
of the surface 12 C abundance. In other words, the LOPSE 
is represented, to first order, by a straight line in the 
[C+N+0]-[ 23 Na] plane with a slope of 3 (thick solid line 
in Fig. pi); 

2. the LOPSE is insensitive, to first order, to / [that pa- 
rameter does not appear in Eq. (|19|)1. In other words, pri- 
mary 23 Na abundance prediction as a function of primary 
12 C abundance is insensitive to the dredge-up efficiency; 

3. the value of 23 X S predicted by the LOPSE is linearly 
proportional to d, as expected since that parameter rep- 
resents the 'destruction' factor of 23 Na by p-capture; 

4. the value of 23 X S is also linearly proportional to the 
parameter ft at given surface 12 C abundance. Indeed, if 
the dilution of the intershcll 14 N into the pulse is low (i.e. 
the dilution factor h is high) then the abundance of 22 Ne 
emerging from the pulse is high, leading to higher 23 Na 
production during the next interpulse. Thus, high ft fac- 
tors imply high 23 X S abundances at a given 12 C abun- 
dance, as expressed by Eq. (]l9|); 
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Fig. 14. Same as Fig. ||, but a for 12 X D = 0.1, 0.2, 0.3, 0.4 
and 0.5 at fixed /, h and d, b for / = 0.001, 0.005 and 0.01 
at fixed 12 Ad, h and d, c for h — 0.2, 0.4 and 0.6 at fixed 
12 Xd, / and d, and d for d = 1.0,0.5 and 0.25 at fixed 
12 Xrj, / and h. All tracks correspond to Z = 0.001 initial 
composition. 



5. the most influential parameter on the LOPSE is 12 Xrj, 
which appears to the second power in Eq. (U9). This re- 
sults from the fact that the surface 12 C abundance 12 X S 
increases linearly pulse after pulse, the 12 X S increment 
being proportional to 12 Ad, while 23 X S increases to the 
third power of the pulse number. The relation is such that 



23 X S increases with decreasing 12 Xu at a given surface 12 C 
abundancen 

Non-zero metallicity stars 

The evolution of 23 X S at Z =4 is dictated by Eqs. (|l6|), 
and shown in thin solid lines in Fig. [l| for several metal- 
licities (the initial abundances are scaled from solar to the 
required metallicity as represented by the dotted line in 
Fig. |l3|). The surface 23 Na abundance stays constant until 
the surface C+N+O abundance reaches the LOPSE. This 
results from the fact that the abundance of primary 23 Na 
produced during the first pulses (i.e. when the C+N+O 
abundance is still far below the value required by the 
LOPSE) is much lower than the initial 23 X S . No variation 
of the surface 23 Na abundance is thus expected. Only when 
the surface C+N+O abundance approaches the LOPSE 
does the increasing contribution of primary 23 Na become 
competitive with 23 X S . From then on, the surface 23 Na 
abundance is mainly provided by the primary 23 Na, and 
23 X S increases asymptotically towards the LOPSE. Those 
conclusions confirm the ones obtained from the synthetic 



predictions presented in Sect. 1.1 



Sensitivity to parameters 

Sensitivity to 12 Xrj The sensitivity of 23 X S predictions to 
12 Xrj is illustrated in Fig. 14 1. As already discussed in the 
analysis of zero metallicity stars (see above), 12 Xo has a 
main influence on the location of the LOPSE. As a con- 
sequence, it affects directly the evolution of 23 A S at all 
metallicities. The lines at different 12 Xd values are practi- 
cally parallel to each other, the ones with the lowest 12 Xd 
predicting the highest surface 23 Na abundances. At a given 
[C+N+O] value, the increase in [ 23 Na] on the LOPSE is 
twice the decrease in 12 Xr> expressed in dex [ 23 A S is in- 
versely proportional to the square of 12 Xp in Eq. (|19l) 



Sensitivity to f Figure \\Ap reveals that 23 A S is almost 
insensitive to /. This property is expected from the inde- 
pendence of the LOPSE on that parameter [see Eq. (fOj)]. 

Sensitivity to h The effect of h on 23 X S predictions is il- 
lustrated in Fig. p+lc. The importance of that parameter 
results directly from the linear dependence of 23 A S on h 
on the LOPSE. At a given [C+N+O] value, the increase 
in [ 23 Na] on the LOPSE is equal to that of h expressed 
in dex [ 23 X S is linearly proportional to h in Eq. (|l9|)l . We 
further note that the value of [C+N+O] from which 23 A S 



9 The fact that 23 X S decreases with increasing 12 Xd seems 
against intuition, since 23 Na is synthesized from the inter- 
shell 12 C. What actually happens is that the surface 12 C, and 
thus [C+N+O], increases more rapidly with 12 Xd, shifting the 
tracks in the [C+N+0]-[Na] diagram to the right. 
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begins to increase is independent of h, contrary to the case 
of 12 X D (cf. Figs. 0a and c). 

Sensitivity to / The effect of d on 23 X S is similar to the 
effect of h, both of those parameters influencing linearly 
the LOPSE. Its effect is illustrated in Fig. 117 
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